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by 
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ABSTRACT 


This thesis is a theoretical and experimental study of a 
Scheme for conve,ting energy from a streaming fluid through an 
elastic plate to an electric circuit. The model is a quasi-one- 
dimensional motion of the elastic plate under the influence of 
the supersonic fluid and a transverse magnetic field. 


An analysis of the model leads to the equations of motion of 
the plate in the unrestrained case. By assuming traveling wave 
solutíons, the díspersion relation, resulting from the coupled 
equations of motíon, is solved. The boundary conditions are dis- 
cussed and the resulting eigenvalue problem is stated. By the 
method of coupling modes, the effect of the coupling on the wave 
numbers is determined. 


In the experiment, an elastic media is tested for response 
with electromagnetic and mechanical excitation. The supersonic 
flow is attained with a set of convergent-divergent nozzles 
feeding a test section. The combined plate end flow experiment 
is run and results discussed. 


The conclusions are: (1) the theory predicts an instability 
in the natural modes; (2) the plate can be excited in the second 
modes and the frequency response íncreases with tension. (3) The 
effects of the instability in the combined experiment are masked 
by forced excitation probably due to the turbulent boundary layer. 


Thesis Supervisor: James R. Melcher 
Title: Associate Professor of Electrical Engineering 
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CHAPTER 1 


Introduction 


Energy conversion has historically been accomplished by rotat- 
ing machinery, making use of small gap lengths and conduction in good 
conductors. Recent investigations have lead to novel static and 
dynamic converters. The feasibility of generating d.c. power from 
the interaction of a moving conducting fluid and an electromagnetic 


a,» Several schemes for A.C. 


system has already been demonstrate 
power generation have been presented, however, these are restricted 
by the requlrement of large magnetic flelds and a good conductine 
fluid. An alternative method for converting energy is shown in 
Figure 1.1. It is a device for producing A.C. power by abstracting 


energy from the fluid flow. 


The devíce consists of an elastic plate in a channel with a 
fluid flow on both sides. The plate is constructed of highly per- 
meable material and is free to move in the transverse direction. 

The channel, also of highly permeable material, contains the fluid 
and copper conductors. One set of conductors is excited by a con- 
stant current (D.C.) source, causing a transverse magnetic field. 
The interaction between the fluid and the plate causes transverse 
motion of the plate, similar to the '"flutter" motion in aero-elasti- 
city. This motion provides a variable coupling of the first cir- 
cuit with the second magnetic circuit, producing A.C. power, as 


indicated in the lumped parameter model shown in Figure 1.2. 


Analytically the problem of the thin airfoil performing small 


lateral oscillations in a uniform stream of incompressible fluid 


Jee 
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FIGURE 1.2 


was for many years the heart of all flutter prediction. The com- 
plete solution was first published in the United States by Theo- 

This approach has been further developed by many workers‏ و وروت 
and is presented most comprehensively in the text, Aeroelasticity,‏ 


by Bisplinghoff, وو‎ )2( 


The definition of panel flutter is the self excited oscilla- 
tion of the external skin of a flight vehicle when exposed to an 
airflowfon one side ‘The analysis of panel flutter has been 
pursued in great detail recently, but there is still some disagree- 
ment between theory and experiment. It appears that the theoreti- 
cal considerations to date do not approach the situation described 
here. The tensions on the panel are significantly less than those 
to be applied in this device. The mechanism for flutter is not 
explained physically in the context of a feedback situation. How- 
ever, the effects of aspect ratio, elastic foundations, damping 
and boundary support have been covered in detail. In the appendix 
a comparision of the prediction of flutter theory and the pre- 
dictions of this theory is shown. The aeroelastician in general 
has confined his analysis to that of "infinite" flow, whereas in 
this scheme the confining effect of the rigid walls will be impor- 
tant. In dense fluids, there has also been investigation of flow 
induced vibrations of fla 
follow: 
recent; 


6 
we) The latter is concerned with a statistical repre- 


elastic. 
sentation of the turbulent boundary layer as a forcing function. 
Here arain, the models have infinite geometry and, therefore, are 
not directly applicable. Lastly there are several "handbook" solu- 


tions available for determining natural modes of vibration, for 
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(7). 


example, see Den Hertog 


Melcher ۰ has studied field coupled surface waves for both 
electrohydrodynamic and magnetohydrodynamic systems. Ketterer(?? 
studying under him, investigated the interactions between a fluid 
stream and a fixed structure capable of supporting traveling 
waves By theory and experiment, it was shown that energy could 


be extracted from the flow in the form of oscillations. 


In this study there is contact between the fluid and the 
elastic structure which was not present in the experiment conducted 
by Ketterer. The effect of this contact, a boundary layer which 
could be turbulent, may be considerable in comparison t^ the 


effect of an instability in the natural modes. 


There are several methods of exciting the plate available. 
First, it could be located transversely in the flow path, simi- 
lar to holding a piece of grass between the thumbs and blowing. 
However, the length of interaction would be limited and therefore, 
the amount of energy taken from the flow would be limíted. 
Secondly, the plate could be fixed at one end only, and the exci- 
tation would be dependent on the trailing vorticies. This would 
be similar to the study reported by efe. Another method 
that has been suggested is to place the plate in a tube closed at 
the downstream end and put the combination in the supersonic flow. 
This would provide greater excitation since pressure variations 
considerably above stagnation pressure have been observed in such 
a configuration. The overall efficiency of this scheme is ques- 


tionable, however. 


The method chosen was considered both theoretically and experi- 
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mentally the simplest to study, with the added feature that the 
significant interaction length might oroduce a reasonably effi- 
cient idevice. In Chapter TI, the device is studied by assuming 
a quasi-one dimensional motion of the plate. The resulting dis- 
nersion relationshin is solved in .hapter I[l and the imolica- 
tions discussed. In cha>ster IV the longitudinal boundary condi- 
tions are treated. The exveriment and results are described in 
Cha»ter /. 
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CHAPTER TI 


ANALYSIS OF DEVICE 





The device will be modeled by assuming a quasi-one dimen- 
sional motion of the olate. The model, as seen in — له‎ 
from a top view, is shown in Figure 2.1. The flow will be assuned 
inviscid and Ancompressible in this analysis since the objective 
is a device that can be tested in a wind-tunnel. The »roposed 
orientation vill also allow neglecting the effects of gravity. | 
The fluid 1 the olate and thé magnetic effects wili be treated 
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for an inviscid fluid and the equation for continuity of mass is 


ወር + ° = 2 
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for compressible flow We will assume small perturbations and 
that the waves are long compared to the significant dimensions, 


that is: 
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ANALYTICAL MODEL 


FIGURE 2.1 
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med quantities the "perturbations. -inserting these into equations 


(2 1) and (2.2) and retaining "^ first order terms: 
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Initially, the thermodynamic process will be considered adia- 
batic, which 8115. gelirine the changes in pressure and density 


by a constant. ۱ — — — | 
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3p = on | 








2.2 Magnetic Equations 
| 


We assume a perfectly magnetic material so that the coil will | 
| 


| 
| 
! 
7 A section of the magnetic circuit is shown in Figure 2.2. 
| induce a magnetic field in the circuit, with the H field existirg 
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in the gap only. From Ampere'S circuital law we get: 
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Define: Hy = + as the steady state H field. Therefore, the 


magnetic fields with reference to figure (2.1) will be: 
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where the superscripts refer to the resnective sides of the plate. 
There are several methods available for computing the force on 
the plate. The most direct method is to use the Maxwell Stress 
Tensor described in detail in A pta. The force lensity 

in the y direction acting on the plate is the sum of the stresses 
acting in that direction on each of the faces of the plate. The 


stress tensor: 
Ih 
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Therefore, the force is proportional to the displacement and will 


have a destabilizing effect on the plate. 
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The plate A" a denfity p and a modulus of elasticity E. 
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Now, summing the moments about the lower left hand corner: 
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The three previous developments can be combined to give the equa- 
tions of motion. 


2.4 The Equations of Motion 


Equation (2.4) ean be rearranged and integrated directly. 
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The transverse fluid velocity must be zero at the channel wall 
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If one assume8 no cavitation the transverse fluid velocity must 


equal the plate velocity at the plate. 
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the above two equations maybe added: 
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It*will be convenient to have the equations in nondimensional 


form before further manipulations: 
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We have normalized the oarameters to a dynamic pressure 
E p v^) and then multiplied by the density ratio. This will 
simplify the solutions as well as provide an indication of the 
$trength of the various forces or pressures compared to the fluid 
dynamic pressure. 2 
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The above form the equations of motion of the system in nondimen- 


sional form 


To solve these equations as they stand for specified forcing 
function would be difficult indeed. This approach has been fol- 
lowed on simplified equations, where the forcing function was a 
statistical representation of the boundary layer excitation. The 
approach we will follow is to look for instabilities in the free 
modes, when the source of oscillations is not a lossy one as it 


is in a turbulent boundary layer. 





CHAPTER III 
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| THE DISPERSION EQUATION 


The dispersion equation, or the relationship between w and k, 
will give an indication of the stability of the plate in the free 
modes. We will assume a delta function in time and space as the 
excigation, so that the equation of motion will have no driving 
function. In Brder to determine whether the structure can suv- 
port traveling waves, we also assume that the solutions will be 
of the form, = 4 ,J (ot -kx) anda. d $ ei (wt-kx), Insert- 


ing these into the equations of motion we get: 


Hee? -ል"ንዳ « ODE = 0‏ نس رد 

(3.1) 
-k (8-9፪*. ነቴ + @2)8 + Ë دم‎ = 0 

(3.2) 


Setting the determinant of the coefficients equal to zero 
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This is the desired dispersion relationship, and its solu- 
tion, o = £(k), and the inverse, k = f(#), will provide insight 
into the stability of the infinite plate. Before proceeding with 
the solution, there are several limiting cases which should be 


investigated. These in effect will separate the two systems. 


1) If we allow the tension, P, on the plate to become 
infinite, the plate becomes,a rigid boundary. 


Thus, if ج ا‎ , we get: 


0 = وت سس + GEI‏ 


or, in dimensional form: 


: 2 _ ¿ap? | 
(w- VK) = Mk (3.4) 


which is the expected dispersion equation for convective acoustic 


waves 


2) Similarly if we allow the density of the fluid to approach 


zero, we find the equation of the plate alone. 
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It is comforting to note that the magnetic force will have a 
destabilízing effect, while the tension and rigidity tend to sta- 


bilize the solution. 


3) Allowing the speed of sound to become infinite will res- 


trict the solution to the subsonic region. In dimensional form: 
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This is the equation for the plate with a dynamic pressure term 


which can be derived from Bernoulli's equation in a linearized 
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The above three checks indicate that the model and the ana- 


lysis are consistant. It also anpears that they are to some 
extent physical. o JP کا‎ TM MAT, መ. 
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The uncoupled equations (3.4) and (3.5) can be solved simply, 


since they are quadratic, and plotted, as in Figure 3.1. For the 
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rest of the analysis, some experimental considerations will 
govern the approach. Supersonic speeds are difficult to attain, 
and in general, the experimental apparatus is designed for a 
specified flow Therefore, we wil) take م۷‎ to be a constant 
equal to twice the speed of sound. The rigidity of the plate is 
also.a fixed quantity for a given experiment and the effect of 
its variation would be roughly equivalent to varying the tension 
on the plate. Since the plate tension is easily changed experi- 
mentally, it will be allowed to vary in the solution of the dis- 


persion equation. 


We now rearrange the dispersion equation as a polynomial in 
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and, using the computer, solve for real values of k. 
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The results are shown in Figures 3.2, 3.3, and 3.4 for the 
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following cases: | ات‎ 


1) Figure 3.3 | M = ,675 (strong tension) 
2) Figure 3,3 1 = ¿545 (medium tension) ` 
3) Figure 3.4 M = ‚314°; ~~ (weak tension) 
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The solutions plotted-ín Figure.3.1 remain as the asympotes 
for these solutions. Retwrning to equations 3.3 and 3.4, we find 


in dimensional form: 
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In the case of strong tension the coupling is weak. However, as 
the tension is decreased as in Figure 3.3 to a point where vr 
the wave speed in the plate, is approximately equal to (V-c), 

the "slow fluid wave, the coupling becomes pronounced. Decreas- 
ing the tension further, decreases the coupling until the solu- 


tions are uncoupled 


The two systems, the fluid and the plate, are in themselves 
stable at at least if the majae tic field is small. The above 


solutions indicate that when coupled together an unstable situa- 
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tion is produced. That is: œ has a negative imaginary part 

and the assumed solution: nan e] (ot- -kx) has the form: y = ፈርን መ 
such that the displacement grows in time, This is the result 

of a feedback mechanism within the coupled system. Let's assume 

a deflection of the plate in the second mode, with a positive 

amplitude over the leading half of the plate. This provides a 

compression wave in the fluid. As this disturbance is swept 

downstream, it arrives over the trailing half of the plate in 

time’ to reinforce the motion of و‎ Ther the motion is trans- 


mitted 


TN‏ ظ 


3 


back upstream through the plate, to arrive at the right time to 
start the procedure over again. This is the case when the velo- 
cities سا‎ and V-c are equal. Since this is a distributed pheno- 
menon over the half length of the plate, the mechanism will work 
for velocities that are approximately equal as well as exact!y 


equal. 


A complete discussion of the dispersion eauation as it per- 
tains to the infinite problem would not add significantly to the 


study of the device proposed. However, the steps will be outlined. 
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Having solved for w = 8 ሐ) the complex values of w obtained 
are substituted and the eguation solved for k. * jk; = f (o, it jo) 
to determine the type of instability involved. This can be done 


(11) 


using the Bers-Briggs criteria, with the help of computer 


technique giving Ss graphical display of both o versus Ds and 
k, versus k, as described by Kills, (1? 
The dísoersion eauations shown in Figures 3.2 - 3.4 have 
eharacteristics which indicate that the basic theoretical model 
is physically meaningful. The coupling between the flow and the 
elastic media occurs at long wave lengths, with the short wave- 
lengths stable and uncoubpled. These are the nronerties of the 
stream structure nreviously investigated in the context of an 
electrohydrodynamic ٠۷١۴ ۵و ذذ‎ P? There it was shown that the 
effect of the longitudinal boundaries was to make the instability 
absolute. That is, the theory supports the conjecture that grow- 


ing oscillations cari be produced. 
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CHAPTER IV | 


BOUNDARY CONDITIONS 


| 
| 
| 
| 


The application of the boundary conditions will be approached 
. by two MTS exact method resuiting in the eigenvalue pro- | 
blem and an approximate method called the coupling of modes | 
approach. A sigple linearized examoie will be covered first fold 
| lowed by a discussion of how to arly the boundary conditions. | 
| This leads to the statement of the eigenvalue problem. Next, the 
coupling of modes is discussed. This will give an indication of 
a more meanings way-of agplytng-the ۳۲ sodfretoRe, and 


| an estimate of the couoling effect. The boundary value problem 
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| is not solved completely. 


41 A Simple Example 


Consider a membrane fixed at the ends. This is similar to 
the elastic plate if "ee یر‎ 2 and I ی‎ are neglected. 


A MGM " NP 
۱ 2 2 
1 Equations of motion: QE. 4 d f 
| At^ dx? 
| 2, Assume Solutions of the form: € =< (A sin k. + B cos k.) m 
0 a 2 


0 « 2 |- + YE T 
3. Dispersion RelationsHis: — w" = b/f k 
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4. Boundary Conditions: Tw" ኝ 


۲ Met) = በ B= 


L 


ds | ጻዴ EL, U= uh 0 kL = D and k = = 


5. Solution: 


1 
۱ ۱ ` 2 
6,» - ٥٥ É (an مزه‎ EO cS mot] 


4.2 2. Method of Characteristics 


In the example the basic steps to the >roblem solution are 
stated. There is a significant change when the fluid is added, 
especially when it is allowed to have a gross translational 
motion The question arises - at what point may one apply boun- 
dary conditions on the flow such that they are meaningful? Were 
the fluid motionless, the boundary conditions must be applied 
at two points in space. This is also true if the velocity of the 
fluid is less than the speed of sound in the fluid. From the dis- 
persion equation, however, the instabilities oceur at supersonic 


speed only. 


¿An order to understand this more clearly, the method ^‏ ہے 
characteristics is useful. | This method uses the general solutions‏ 


of the wave equation for the transient case: 


f = f + (x-ct) + £ - (xt ct) 
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| 
| 
| 
| 
| 


where f+ is a vaye-propagating in the-plus x direction and f- a 
wave in the negative x direction. The SION ona are represented 
by tbe third dimension. of a. graph where the first two are time and 
space. Thus a disturbance can be followed 49- both time'and space. 


x 
| 
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Let : Rég ë and Dea = c 


— u .- ج ت‎ p... 


Assume a disturbance: p = I - a ¿< KK < + Bü 


— 


0 elsewhere. 


for the fluid flow shown in کم .4 د راس‎ 


m — — 


In Figure 4:3 (ከ5 9850258668 ig shown. ` ۱ | 


Lines representin& constant 2 and? can 'be drawn on the 


graph Those originating at.x = + a are\the most interesting 


— — — 


since they show the progress of the disturbance. For example, 


at t = ed the #፦፦ of NA distu mom. cin be determined from: 


ay 


E t 
O 


- ይ፥(5) ቪ - E + ይ.“ 


رش 


/ 
t 


The siene of the lines is of interest since, if the lines of 5 


had a vositive slope, the point **= -L, would never know the dis- 


turbance This is the case when the fria is allowed to have 


motion at a velocity greater than the speed of sound. For a trans- 


lating media, we must redefine: 


4 
wa <- — _ 
- - < 


í = x +(v+c)t and 3 = x+(v-c)t 
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as the characteristic solutions to the convecting wave equation. 
Then (v-c), the slope of the 9 lines, becomes positive for v > c, 


and points upstream will not be effected by events downstream. 


This implies, in the scheme proposed, that the fluid boun- 
dary conditions must be applied at the leading edge of the »late 
in order to dette hine their effect Jon the interaction. The fluid 
will not know what is happening at the trailing edge until it 
gets there. There are two variables on the flow, the pressure 
and velocity, and their values at the voint x = 0, as a function 


of time, form the constraints on the flow. 


4.3 Ihe Eigenvalue Problem 


The dispersion equations was sixth order in œ . The six 
boundary conditions are the two specified for the flow, plus the 
four on the plate. From simple beam theory the latter are deter- 
mined from the tyoe of support-free, pinned, or clamped. If we 
assume the ends are clamped, this implies that that the deflec- 
tion and slope, s, are zero at x = O and x = L. Since v(x,t) 
and p(x,t) are definel as difference quantities (above and below 


the plate), their values are zero at x 20. 


Returning to the equations o£ motion this can be stated sim- 
ply. 

^ - pe u 2-2. ^ 

(u - k) (uMa-P.E.)n (٢ 2 (Ma عا‎ ( v = 0 (&.1) 


K(P.E.)3 «Y (b-k) v 99 (4.2) 
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(B.E.)n t. p-0, | (4.3) 
where P. E. - وج ده دنو‎ 


As in the example, assume a solution of the form: 


6 ۲ 
n(x,t) = Re دا‎ A: e JE x, J8nt ) 


r 
| 


| and rewrite equation (4.1) and (4.2). 


1 


. B: 





۵ = - zur A | (4.4) 


v ke CH. م‎ ۶ 
k(P.E.) 


applying the boundary conditions: 


Di n(0,t) =.0 f Ae و سر‎ (498) 
- ፳5 (5,5)-0- 2 ፡%ቂ = 9 (4.65) 
(3) و‎ < A n= 0 N (4.60) 
Wirt ae -0 هماند‎ eMe wp (4.64) 


(5) -p(90,t) = 0 £ 20 au Á r = 0 (4.6e) 
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(6) v(0,t) = 0 E 2£, (k_) ል =ዐ (4.6፻) 
7 | | n 


Rearranging in matrix form: 


ለ(,፪)= ' 

1 1 1 1 1 1 

k, k, k, ky k, kc 

WU | "LM "LOI. ۱ ak. m AP ei Bel 

E.M ` k ei لو دي چ‎ k a و‎ k حم‎ k_ ed 'ነ; 

-1 2 3 EA 5 6 

7 =..=ጩጪሕ غه‎ — 1 = 1 — - 
P.E (kK, P.E (ko) P.E.(k.) P.E (K) 2-8-0717 P.E. (Kg) 

£4 (k4) £, (ky) £, (k4) f (k,) £; (kç) | f, (kg) 


In theory then the equations ^ (ay) - 0 and D(o,k) = 0. can 
be solved. Assuming a complex w, the six k's are found from 
D(w,k) = 0. These values are substituted in to ^. This procedure 
is continued until ል =6; both it's real and imaginary parts, 

Since it is complex In practice another method, that of weak 


coupling of modes will be used to attack the problem. 
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4.4 Method of Coupling of Modes 


One might ask the question: "Tf two systems are individually 
Stable, why would a weak coupling between the two result in totally 
different motions?" The answer i5 that the basic solutions for 
the two systems are approximately correct, and that the changes can 
be predicted within certain limitations. First, the systems must 
be stable by themselves. [In this problem, this restriction will 
not allow . including the effect of the magnetic stress on the 
plate since it has destabilizing influence. Secondly, the coup- 
ling between the two systems must be weak. Within these restric- 
tions, mutual coupling of modes will allow simplifications based 
on a physical understanding of the system. 


(13) .. A- comp ita Ma 


In theory, the approach is quite simple. 
system is first divided into its component parts. The equations 
of motion for the isolated systems are solved exactly, and the 
solutions are expressed in the "normal modes of the element. 

The complex system is then assumed to be made up of these ele- 
ments weakly coupled. The motions within the system are described 
by the elemental motions which are slightly perturbed by the 


coupling. 


As an example, in the system under consideration the ele- 
mental parts are the fluid and the plate. Assuming the fluid 
motionless and the plate of negligible rigidity (a membrane), 
each will propagate disturbances as waves, both a forward wave 


and a backward wave in each system. As seen from the solution 
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of the dispersion relationship, the coupling will occur between 
the forward membrane wave and the "backward’ fluid wave. There- 
fore, the other two waves can be neclected. The motion of the 
system will be that of these two waves slightly altered by the 
coupling. 


4.5 Elemental Equations in the Coupled Mode Form Žž _ 
1. Fluid Equations 


The fluid equations can be greatly simplified to get an 


approximation to the solution. 





:.» ۸۰ سے‎ Ap(x,t) 
— نت‎ eg (vite? At s. 2 | (4,8a) 
Mo(x.t) ` yv (g. t) 


"x "Po dt (4.8b) 


These coupled equations are most easily solved if a method of 
decoupling them can be found Following, sour, (9) the 
decoupled form is derived from linear combinations of the above 


equations. That is; assume a normal mode form: 
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where the mode amplitudes are normalized so that the square of 
the amplitudes represents energy carried by the mode. 


1 
Po (Vtc) ) 


and adding and subtracting from سن چا‎ (A. Bai will Em Ce 
desired — a form: ۱ ۱ ' | 


Multiplying equation (4 8b) by the fluid impedance ( 


u F ١ ١! - A 
7/8 MET. Q Q Ql e done 
د ا يو‎ p (Vte) — ^^ mr w ۳۹ 
or 
bor pg په لس وا‎ (34:07) 
= Vte Qt = په‎ 


Assume that the time dependence will be of the form er 


jut 


then: b + (x,t) = a T (x) e A a*- (x) @ Jat 


Now if v(x,t) = Re oe 


then a + (x) -3 9ا 2ری‎ + — P(x)) 
ga 


and equations (4.8a) and (4.8b) become 


) 2 + رم و مل د رر‎ 6 (4.9) 
b 1 4 ( A 
( دا‎ jk)a-(x) = 0 da ji9a* (x) = 


(4.9b) 
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2. The P'ate Equations 


Again simplifying assumptions must be made. 
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Neglect- 


ing the rigidity and the external forces on the plate, the nor- 


mal mode form of the equations for this element can be derived 


in an analogous manner. 


-Av (x.t) Ti ። — where v(x,t) ¿2 


8... VS (x,t af 


then 


ም፦ 


Bt (x,t) = 2 (p 2(VGx,t) + Š S(x,t)) 


and A 

w Ca 

ጫኗ (b + Crt e (HS (bt (x,t)) 
there fore, 
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where 
0 


ty = سه‎ 
5 (P/p)* 


4:6 Coupling Between the Two Elements 


In section B, we derived the normal mode form of the element 
equations In effect the single variables in equations (4.8a) and 
(4 8b) have been replaced by a combination of the two variables. 
To put the previous equations in coupled mode form, it is nec- 
essary to make linear combinations of the elemental forms. This 
is done in such a manner that when the coupling is reduced to 
zero, the equations reduce to those derived for each element 
separately. 
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4 a? de | + 
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In order to make the problem more tractable we must make use of 
our previous results. The first is that coupling between the two 
systems will occur when the two wave members, ky and ko; are 
approximately equal. Since k, - — and FI and in Chapter 1 
we saw that the coupling occured when ፍ፡ ተ= ሚወ we can there- 
fore neglect the fast fluid wave (V*c) and the backward plate 


wave « E 


Secondly, we must determine if, in fact, the coupling between 
the two systems is weak. Returning, to the coupled dispersion 
relationship and inserting our previous assumption that the 


rigidity and the magnetic field are zero, we have: 
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or in dimensional form and rearranged. 
fe 2 | Lei, 2 _ ۸ (or ky? 
0 2 m. © 32 0 


where‘ the term on the right represents the coupling between the 
two systems: Now let the perturbation of the elemental wave 


numbers be small: 


Che 5. 
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L The coupling is. proportional to ratio of the densities 
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Therefore, Coo = jk, and equations 4.11 reduce to: 
"es i ik + ° * ( 
مس ارا‎ ብንን | | (4. 128) 


Ax - 53252 + Jka (4,125) 


1 e 


١ 
| qi i á | 


lhese eguations can be solved by assuming solutions of the form 
e * ¿and setting the determinant of the coefficients equal to 


Zero. 
6 ሐሑ 
F kk) 1*3 
TERA - ጫ“*ቃቻሽ ጥገ 12), (4.13) 
where k, = چا‎ and ko = 8 the undisturbed elemental wave 
numbers. S 


Since they are approximately equal: 
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4.7Summary 


The eigenvalue problem has been stated, but was not solved 
because the method of coupling of modes gave better insight into 
this problem. We applied coupling of modes to a simplified pro- 
blem and determined the perturbation of the wave nembers based 
on ‘the effect of the coupling. In order to properly apply the 
boundary conditions in this case, the approach of coupling of 
modes could be followed. But it is not clear that the neglected 
waves can in fact, be neglected. Since at long wave lengths, the 
boundaries will have a definite effect on both modes, this must 
be accounted for in the theory. The results of Chapter III and 
IV do provide a basis for experimentation. If we assume that 
the wave numbers are only slightly perturbed; we can compute kL 
from the physical length, k, from equation 4.13 and w and 7 
from the dispersion relationship. Although these values are not 
exactly self-consistent, they should be approximately ‘so. 
Experimental techniques to achieve these values are described in 


the next section. 
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CHAPTER V 
EXPERIMENTAL PROCEDURE 


In order to experimentally test the predictions of the 
theory developed, it is necessary first to insure that the 
"ingredients" required are available. The initial attempt at 
assembling a device and making test runs was not successful. 

The step by step approach will be described in detail. This con- 
sisted of several séparate experiments. The first was testing 

of the plate structure to be used. The second was investigation 
of a means of attaining supersonic flow. Next, the two were com- 


bined to test the flutter theory. 
5.1 Testing of ‘the {Plate 

There ato 77 performed on the plate along. The 
first was an investigation of various types of material and vari- 
ous end conditions Thertest set up is shown schematically in 
Figure 5.1. Due to the m@enetic coupling between the two cir- 
cuits, the results of this test were qualitative rather than 
quantitative. The results are shown in table 5.1. The end 
conditions have a noticeable effect on the response. The damp- 
ing or losses #2Ycondition but tightly clamped appeared to nulify 


any response. 
The frequency did increase with tension which indicates that 


the use of membrane analysis’ ‘in the coupling of modes section is 


not totally in error. The frequency of response was reasonably 
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TABLES 
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close to that predicted. (See example 5.1) 


-» 


۱ Ti 2 
EXAMPLE 5.1 اک سو‎ n 
Cer و‎ ú 
where -— 
D = EL, the rigidity 
u = mass per. unit area | 0 
M = longitudinal harmonic mode * 
° = Length 
5 1 ! | A 1 1 
2 650 በፍር | 71 
£, = 104 eps 


e. ድ‏ 1 ب 


| 


The second experiment on the plate above was to determine 
if in fact the second harmonic frequency could be excited when the 
media had the perdominate properties.of a membrane. Equipment 
similar to that shown in Bette 5.1, was used for this part of 
the experimentation fo avoid the magnetie coupling, one of the 
electroragnets was replaced by a speaker for machete ett excita- 
tion with weak coupling. The plate thickness was reduced to .001 
inch, to more closely.approximate a membrane. ‘The equipment is 


shown pictorially in Figure 5. 


In thts configuration the plate frequencies were not detec- 
table The membrane frequencies pave excellent response and are 


tabulated in Table 5.3 The variation between predicted and 
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actual frequencies can be explained. First, the actual para- 
meters for the material are difficult to determine exactly. 
Secondly, at the higher values of tension, there is a possibility 


that. the grips or end clamps slipped slightly. 


The maín conclusion to be drawn here is that the second 
harmonic can definitely be observed. The quality of the response 
decrease with increasing frequency since the higher modes are 
more difficult to exert. The frequency expected when operating 
the wind tunnel could be reasonably predicted from this type of 


experiment 


The responses at 355 and 720 cps are probably due to weak 
coupling with the speaker or other portion of the apparatus, since 


they are independent of tension. 


5.9 Experiment to Obtain Supersonic Flow. 





The physics of supersonic flow is a very complicated sub- 
ject and is extensively covered in the literature, notably in 
reference 15. Only those princíples which are germane tothis 
experiment will be touched on. In general, supersonic flow can 
not be attained directly from an air compressor. It is necessary 
to modify the compressor output by passing it through a set of 
nozzles. For this experiment the nozzles were of the convergent- 
divergent type. The design procedure is to attain Mach 1 flow at 
the throat of the test section and then expand to the velocity 


desired. The calculations for this are shown in example 5.2, 
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based on the set of nozzles shown in Figure 5.4 


EXAMPLE 5.2 (Reference 15, Table BITT) 


Inlet Conditions: 


vsia T = 70°F‏ 125 = م 
o 0‏ 


a. Static Conditions: 


P 
0 mn 
p ao A. OS E 
5 T, £t? 
c - 49.1 T = 1158 - 
0 0 sec. 


b. Throat Conditions: 


P*/P = 0.528 


0.823 - ری 

p*/p, * 0.634 

ር*/ር = 0. 913 
ዕ 


ር Test section conditions for Hach No. = 2.0 


2 
= x = = 
Mk = 1.633 A/A 1.687 = 6 , ዐርዉ 
5/8۵ = 0.1278 V — 1733 ft/sec A 4 2 
throat- 3.57cm 
P_ = 16.1 osia 
P1/ Po - 0.2300 L 


TT, = 0 555 


d. Mass Flow 


% يپ 


۵ x60 = 245 cuft/wmin. 
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Secondly, the flow in the test section will be greatly 
effected by the build up of the boundary !ayer. This has the 
effect of a convergent nozzle, which tends to reduce to flow 
velocity toward Mach 1. In the case of Mach 2 flow, the velo- 
city will be reduced to Mach 1 in a length of 31 diameters. In 
thís experiment, the diameter or depth of the channel will be 
taken as one centimeter to allow sufficient magnetic coupling. 
Therefore, the length must be less than 31 centimeters; an ini- 
tial length of 20 centimeters was used. (From this length, the 
thickness of the plate for a desired natural frequency can be 
computed as in example 5.1). Based on these calculations, the 


experiment shown schematically in Figure 5.3 a»! was eonstructed. 


The pressure readings were taken with both the plate in 
position and removed. The readings were taken with a guage and 
compared with a monometer for accuracy. The readings obtained 


are shown in Figure 5.4. 


These results are at a variance with those predicted in 
example 5.2, but they do indicate supersonic flow. The variance 
is due to several factors. First, the ideal area ratios are 
significantly diminished by the boundary layer to an "effective" 
area ratio. Secondly, the capacity of the air ejectors was 
insufficient for running the experiment as desired. Also, the 
pressure measurements were made at the wall of the channel, pro- 
bably in the boundary layer, and not in *he center of the flow. 
The latter problem is difficult to circumvent, and in previous 


experiments (reference 15, page 135) good agreement between actval 
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wall pressure and theoretical mid-stream pressures has been 
attained. Using the lower pressures recorded, the speed of the 
flow can be approximated. Based on the readings at 35 psia inlet, 
the flow is Mach 1.47. 


5.3 An Experiment to Test the Flutter Theory 


We have determined that the ingredients of the theory are 
attainable. The flow can be „ade Supersonic, the plate can be 
excited in the second mode, and the frequency of this mode varies 
with the tension. Thus the results of Chapter IV indicating 
that the perturbation on the wave nımber k can be estimated from 
the coupling seem reasonable. The rigidity tem Dk? can be 
neglected compared to the membrane term M at the frequencies in 
question. This is especially true if the very thin (.001") plate 
is used. The first experiments with the combined ingredients 
used a .015 inch plate. The equipment consisted of the nozzle, 

a plate, a magnetic detection device and an oscilloscope. The 
presentation on the scope was confused, showing standing waves 
and traveling waves and these responses were at several fre- 
quencies. These were probably excited by the boundary layer. 

The tension was varied and the frequencies of the response 
increased with increasing tension. By varying the pressure the 
response would set in at a certain pressure for a given tension - 
quantitatively a higher tension required a higher pressure. At 

a given pressure, providing it was above the threshold, the vary- 


ing of the tension was limited mechanically. That is, the tension 
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could be increased until the plate failed in tension, and if the 
tension was decreased the plate would fail by excessive dis- 
placement. The latter occured at low values of tension. It was 
felt that nothing of substance to support the theory could be 
gained from further experimentation with this set up. "Flutter" 
did oceur and this did provide a variable coupling between the 


two magnetic circuits and there was an output signal. 


However, this flutter could not be attributed exclusively 
to the mechanism proposed in the theory. It was probably caused 
by excitation from the boundary layer. There was a component in 


the vicinity of the frequency expected as computed below: 


(1) length of membrane: k -tr = „628 TE 
(M -1)kv, 
(2) from dispersion equation: ፻ — — — 1080 eps 


Q, 7 495 ግ 


3 
q= ሜጫ e .25 مې‎ 


(3) from coupling of modes: 2 2 
4d” (v-c) 


In order to prevent the end connections from interferring 
with the flow they were removed from the test section as shown 


in the figure. Therefore, the flow is in transition over a short 
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section at both ends of the plate. This could also provide an 
excitation which was not predicted and would be difficult to 
analyze. 


Therefore, another experiment was designed. Due to a time 
problem, it may not be completed before the due date of this 
thesis. The changes include using a .001 inch plate to more 
closely approximate membrane performance, refined ‘end connectors 
to insure a good response and no flow disturbance by the tension 
mechanism. The results of this phase of the experimentation will 
be reported in Appendix B. 
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CHAPTER VI 


CONCLUS TONS AND RECOMMENDATIONS 


6.1 Conclusions 


The analysis of the model does predict the desired pheno- 
menon. The interaction between the flow and the elastic media 
at long wavelengths will theoretically produce growing oscilla- 
tions. The effects of the longitudinal boundary conditions can 
be approached by the coupling of modes. In the simplified pro- 
blem considered in Chapter IV, the effect of the coupling on the 
wave number was shown. This approach could be expanded to show 


the influence of the boundaries on the wave numbers. 


The experimental portion shows that the separate systems 
can be achieved. The plate can be excited in the second mode. 
The natural frequencies of the plate increase with tension and 
the tension required by the theory is physically realizable. A 
supersonic flow over the plate can also be attained in the labora- 
tory. With the equipment used, the maximum velocity was 1.47 
times the speed of sound. The major source of difficulty. 
g@countered in combining the systems was the excitation provided 
by the boundary layer. This wasrot included in the model or in 
the analysis. The excitation gave forced responses, mainly in 


the first mode, although there was a second mode in the response. 


There are several other considerations which have not been 
covered. With respect to the elastic media, there will be eddy 
currents induced when it oscillates in the magnetic field. 
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These represent a loss in the system which could be significant. 
The Curie temperature will impose a thermal limitation, restrict- 
ing efficient operation of the device to relatively low tempera- 
tures. From the stand point of strength, it ls fortunate that 
ferrous materials, which are required magnetically, also show 


good fatigue properties. 


The most severe limitation on the device efficiency will be 
the losses incurred by the fluid in transiting the nozzle and 
interaction area. It was noted in the experimental section that 
the flow would be supersonic only for a short distance. For 
instance, the Mach 1.47 flow in the experiment would be "choked" 
to Mach 1.0 at the end of the test section (six inches), due to 
boundary layer build up. At the same time the isentropic stag- 
nation pressure would be reduced to less than one half its 
initial values. This represents a considerable loss in efficiency. 
A more refined design would improve on this situation, but there 
will be a definte limiting efficiency attainable. 

The subject of supersonic flow in short ducts is somewhat 
complex and experimental data has not been eg وولو‎ ۷ 
In general, the boundary layer thickness increases with length. 
But in some cases, negative friction coefficients have been 
attained, indicating a decreasing boundary layer. The effects 
of the oscillating plate on the boundary layer have also been 
studied but there appears to be no significant change in boundary 


layer thickness. ( 





= " ^39 S - 








mas 1 1/ጥ8,/-4. ሖ6 


REN NE 2l 
Pon riw کروی‎ መል سسجت کت تو وج‎ ከ 


“urn نس‎ የ18፳23:6፤83 52. 85፻ነጣ8 él 3 caltuzeqo አህ de 


JAU anarei rd اد‎ , ብብዛፒፅ፡፡ bo ambos kuaza 
— · pas ets e ad, ALA er 
299998 5፡6)39፤ ለ5 
ንን ده د‎ = 
ہہ کا‎ — 
رت بس مہم ہہ اس‎ ey 
ER E + 
Gh ates Dass ee (መቹ 
ወ تاس جن‎ Way tele ሆን መ 
id Serum) سم‎ u> ومست‎ ግዜ "A dd S 
ri Lë san gn ggf وج‎ estes w Duss ae «ni fun 


dem wm. cose) vitÀe init مده موه‎ e ei IPY 


ai Û salan ar‏ تہ 6۱ تتح ترود تن 
Pas‏ هماما سے سا Sudhi ajal tee‏ 

BEE Ein) nich ጉና] ዛክፍርርባ ሀ. ርመ 2 
ud end wre وو دس‎ vr aan سید‎ mae ni 4 
alls a seyal دوعس مادنا ۸ اتکی مسب‎ 
Gei saln srat wml yya Lage, Act Ae 
xL A — 0 
pum 































መሼ 





62 


The device itself may not be the most efficient way of 
obtaining energy from the flow. However, it is possible that 
to employ it in tandem with 2 ና. conventional energy con- 
verter would improve overall efficiency. The important conclu- 
sion here is that a mechanism is available for transfer of energy 
between the flow and the media. The mechanism itself is not a 
"lossy" one ~ the inefficiencies are incurred in the scheme 
for using the mechanism attempted here. In other words, this 
paper has been limited in scope to a particular device, but 
further research would lead to other methods of using this con- 


cept for energy transfer. 


6.2 Recommendations 


The study of the longitudinal boundary conditions in this 
scheme could be expanded using the coupling of modes approach. 
The coupling coefficients considered here pertained only to the 
coupling between two waves, one in each system. The coupling 
between the other two waves, both with the other wave in its own 
system and the other wave in the other system were neglected. 

The latter still appears reasonable. But the boundary conditions 
will provide coupling between the two waves in one system. We 
neglected this here because at the wavelonsths in question, the 
rigidity term could be neglected and the plate considered as a 
membrane. In the general case this would not be so, and the area 
of longitudinal boundary conditions in coupled systems has many 


possibilities for study. 
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The combined experiment showed the ovwewstanding effect of 
the boundary layer, essentially msking the desired response 
if present. The boundary layer could be reduced either by 
thermal methods or by expanding the flow in the depth dimension. 
If this were dene, the correlation of experiment and theory 
might be attained and the mechanism studied. 


(9 ) 


electrically to study this mechanism.  Fraize 


used a fluid ¿stream anda spring coupled 
(17) studied a 


device in which a fin , situated in a wind tunnel (sub sonici), 


Ketterer 


was excited mechanically to produce traveling waves. Signifi- 
cant reductions in drag were recorded. Replacing the mechanica: 
excitation by magnetic traction along the length, using the coup- 
ling mechanism investigated here, provides the possibility for 


pumping of or propulsion through a fluid surrounding the fin 


This would be the inverse of the problem studied here, but 
would not have the boundary layer problem inherent in super 


sonic flow. 
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Appendix A 


Ae LO A AR — o MS AMO AM eee کہ‎ 


Flutter is^defined as the self-excited oscillations of a 


(18) This definition leaves some 


airfoil in supersonic flow. 
doubt as to the physics involved, but sínce it is restricted to 
supersonic flow it seems conceivable that the mechanism involved 


is similar to the feedback situation described in Chapter III. 


The problem of flutter is stated as a partial differential 
equation, where structural, inertíal, and also dynamic operators 
are combined to determine the equation of motion for the panel. 
Two approaches are followed from here dependíng on the geometry 
involved. One is to apply the boundary conditions to the plate 
only. Since these are two at each end, a "fourth order" eigen- 
value problem results. 


The second approach is to assume traveling wave solutions, 


pro and assuming real values for k, compute complex values 


(18) has used both approaches and compared the 


11 .ل 2ه 
results. tehen recommends which approach to use based on the‏ 


geometry of the problem consídered. 


In terms of our parameters, the equation of motion is: 
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where the term on the right is the ‘aerodynamic pressure. For 
small M, this reduces to 


and the equation is the same as equation (2.10) except for the 
structured damping term E In our case, the aspect ratio is 
zero and the first approach is recommended. That is: the bound- 
ary conditions on the plate are applied and a solution of the 
form 
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is assumed. This forms the eigenvalue problem, and the soluttons 


for c, , 0, are shown in reference 18. Many parameters have been 


introduced and this leads to lengthy computations to compare the 
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the results of the two theories. For the tension, rigidity, and 
plate thickness employed ín the device, the solutions predict Oo, = 100 
and o. = 1000. The approximation is necessary because the 
plotted data must be extropolated to include the values of ten- 
sion specified. These results show correiation within an order 
“of magnitude of those predicted in Chapter V. ( a : 
7 ). Most of flutter investigation has been restricted 
to determining the onset of flutter and methods for avoiding وناب‎ 
occurance. Again using our parameters in figure 10 of reference 
18, flutter is predicted to occur when the tension corresponds to 
making s" + = v-c. This is surprisingly consistent with the pre- 
vious results. The predicted flutter frequency under these condi- 
tions is on the order of 200, where UN equals the natural plate 
frequency neglecting tension. The previous results predict 
Qe = 10o . 

There is considerable agreement between the two results, 


in spite of the difference in boundary conditions. 
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Appendix B 


The previously described experiment was constructed.  How- 
ever, the improved end conditions to lessen the flow disturbance 
were mechanically unweildly. They did not provide the firm con- 
nection desired. Therefore, no further results of significance 


can be reported. 
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